Analysis of Bell inequality violation in superconducting phase qubits 
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We analyze conditions for violation of the Bell inequality in the Clauser-Horne-Shimony-Holt 
form, focusing on the Josephson phase qubits. We start the analysis with maximum violation in the 
ideal case, and then take into account the effects of the local measurement errors and decoherence. 
A special attention is paid to configurations of the qubit measurement directions in the pseudospin 
space lying within either horizontal or vertical planes; these configurations are optimal in certain 
cases. Besides local measurement errors and decoherence, we also discuss the effect of measurement 
crosstalk, which affects both the classical inequality and the quantum result. In particular, we 
propose a version of the Bell inequality which is insensitive to the crosstalk. 
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I. INTRODUCTION 

In 1935 Einstein, Podolsky, and Rosen (EPR) have 
shown in the classical paper^ that quantum mechanics 
contradicts the natural assumption (the "local realism" ) 
that a measurement of one of two spatially separated ob- 
jects does not affect the other one. This "spooky action 
at a distance" - known as an entanglement - is now rec- 
ognized as a major resource in the field of quantum infor- 
mation and quantum computingi^ The paradox has led 
EPR to conclude that quantum mechanics is an incom- 
plete description of physical reality, thus implying that 
some local hidden variables are needed. 

The EPR paradox remained at the level of semi- 
philosophical discussions until 1964, when John Bell con- 
tributed an inequality for results of a spin-correlation 
experiment,^ which should hold for any theory involving 
local hidden variables, but is violated by quantum me- 
chanics. Inspired by Bell's idea, in 1969, Clauser, Horne, 
Shimony, and Holt^ (CHSH) proposed a version of the 
Bell inequality (BI, the generic name for a family of in- 
equalities), which made the experimental testing of local 
hidden-variable theories possible. The main advantage 
of the CHSH inequality in comparison with the original 
BI is that it does not rely on an experimentally unre- 
alistic assumption of a perfect anticorrelation between 
the measurement results when two spin- 1/2 particles (in 
the spin-0 state) are measured along the same direction. 
Many interesting experiments^'^'^ i^'^d° have been done 
since then. The results of these experiments clearly show 
a violation of the Bell inequalities, in accordance with 
quantum mechanical predictions. 

The BI violation has been mostly demonstrated in 
the experiments with photonsj^s^i^i it has been also 
shown in the experiments with ions in traps^ and with 
an atom-photon system;^ a Bell-like inequality viola- 
tion has been also demonstrated in an experiment with 
single neutronsJ^ Experimental violation of the BI in 
solid-state qubits would be an important step towards 
practical quantum information processing by solid-state 
devices i^^d^d3,i4,i5 Experiments on observation of the 
BI violation in Josephson phase qubits are currently 



underwayji^ii^ Theoretical study related to the BI viola- 
tion in solid-state systems has also attracted significant 
attention in recent years -i^i^i^i^i^i^iiS^ 

In this paper we discuss the Bell inequality (in the 
CHSH form) for solid-state systems, focusing on experi- 
ments with superconducting phase qubits. We study ef- 
fects of various factors detrimental for observation of the 
BI violation, including local measurement errors, deco- 
herence, and interaction between qubits (crosstalk), and 
analyze optimal conditions in presence of these nonideal- 
ities. 

The paper is organized as follows. In Sec. [Til we review 
the CHSH type of the BI, and also discuss tomography- 
type measurements using qubit rotations. In Sec. Illll 
we consider the ideal case and describe all situations for 
which the BI is violated maximally. Sections IIVIIVII are 
devoted to the effects of various nonidealities on the ob- 
servation of the BI violation. In Sec. IIVI we discuss the 
effect of local measurement errors, using a more general 
error model than in the previous approaches . 
Analytical results for maximally entangled states and nu- 
merical results for general two-qubit states are presented. 
In Sec. |V] we consider the effect of local decoherence of 
the qubits. In Sec. IVII we discuss measurement crosstalk, 
which affects both the BI (since crosstalk is a classical 
mechanism of communication between qubits) and the 
quantum result. We also propose a version of the BI, 
which is not affected by the crosstalk. Section |VlT] pro- 
vides concluding remarks. 



II. PRELIMINARIES 

A. CHSH inequality 

We begin with a brief review of the CHSH 
inequality, ''■~>^ a type of the Bell inequality usually 
used in experiments. Let us consider a pair of two- 
level systems (qubits) a and b. Assuming that a real- 
istic (classical) theory based on local observables^ holds 
and there is no communication between the qubits (i.e., 
no crosstalk), the two-qubit measurement results should 
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satisiy the CHSH inequalitjiii^ 

-2<5'<2, (1) 

where 

S = E{a, h) ~ E{d, b') + E{d', h) + E{d', h'). (2) 

Here a and a' {h and b') are the unit radius- vectors on the 
Bloch sphere along the measurement axes for qubit a (b) 
and E(a, b) is the correlator of the measurement results: 

E{d,b) ^p++{d,b) +p — (a, 6) -p+_(a,5) -p^+{d,b), 

(3) 

where Pij{d, b) (i,j = ±) is the joint probability of mea- 
surement results i and j for qubits a and 6, respectively. 

The sum of the probabilities in Eq. ^ equals one. This 
can be used to recast Eq. ^ in the form 



S = 4T + 2, 



(4) 



where 



T = p{d, b) - p{d, b') + p{d', b) + p{d', b') ~ Paid') - pb{b). 

Here p{d,b) = p++{d,b), whereas Pa{d') = p++{d' ,b) + 
p^ — {d' ,b) [or pb{b) = p++{d,h) +p — I- (a, 5)] is the prob- 
ability of the measurement result for qubit a (or 
b) irrespective of the measurement result for the other 
qubit [in the classical theory in absence of communica- 
tion between qubits, Pa(a') is obviously independent of 
the direction b and even independent of the very fact of 
the qubit b measurement; similarly for pi,(b)]. Thus, in- 
stead of the probabilities py , one can equivalently use 
the probabilities p, Pa, and pb, and the inequality ([1]) can 
be recast in an equivalent form^ 



- 1 < T < 0. 



(6) 



Notice that both inequalities ([T|) and (O can have 
somewhat different meanings in different physical sit- 
uations. In particular, the results and "— " may 
correspond to the presence or absence of the detector 
"click" (so called "one-channel measurement"^*^); in this 
case a low- efficiency detector significantly increases the 
chance of the result "— ". Another possibility is the so- 
called "two-channel measurement" , in which the qubit 
states"-!-" and "— " are supposed to produce clicks in dif- 
ferent detectors;^ in this case inefficient detection leads 
to three possible results: "-I-", "—", and "no result". Sig- 
nificant inefhciency of the optical detectors leads to the 
so-called "detector loophole,"— which arises because 
effectively not the whole ensemble of the qubit pairs is 
being measured. This problem is often discussed in terms 
of contrasting the Clauser-Horne^ (CH) and CHSH in- 
terpretations of the inequalities, which differ by consid- 
ering either the whole ensemble or a subensemble of the 
qubit pairs. It is important to mention that in the case 
of the Josephson phase qubits (which formally belongs 



to the class of one-channel measurements) the whole en- 
semble of qubit pairs is being measured, and therefore 
there is no detector loophole (if one avoidsii corrections 
for measurement errors), as well as there is no difference 
between CHSH and CH interpretations. 



B. Tomographic measurements 

In some cases, as for Josephson phase qubits, the mea- 
surement (detector) axis cannot be physically rotated. 
However, instead of the detector rotation, one can rotate 
the qubit state.— Let us show the equivalence of the 
two methods explicitly, using the example of the phase 
qubit and assuming ideal (orthodox) measurement. 

The Hamiltonian of the phase qubit in a microwave 
field in the subspace of the two lowest states in the qubit 
potential well |-0o) a-nd is 

H^^ (;ia;,/2)(|7/;i)(V'i|-|Vo)(V'ol) 

+hn{t) sin{u;t + (/.)(|^o)(V'i| + \i'i){y^o\), (7) 

where h is the Planck constant, LUq is the qubit reso- 
nance frequency, uj is the microwave frequency, fl(t) = 
diQE{t)/h is the time-dependent Rabi frequency, c?io 
is the dipole-moment matrix element, and E(t) is 
the electric-strength amplitude of the microwave field. 
Transforming to the qubit basis (the rotating frame) 
jo) = \ipo), |1) = e"^*|7/;i), neglecting fast oscillating 
terms in the Hamiltonian (the rotating-wave approxima- 
tion) and assuming the resonance condition lo = tOq, we 
arrive at the Hamiltonian 



, hn(t) ^ ^ 

H{t) = ^r?-CT, 



(8) 



where the unit vector ft — (sin 0, — cos 0, 0) lies in the xy 
plane making the angle 4> — n/2 with the x axis, whereas 
f? = {<^x, Cy, az) is the vector of the Pauli matrices. "^^ Here 
and below we associate state |1) (|0)) with the measure- 
ment result "+" ("— ") and with the eigenvalue 1 (—1) of 
az, so that az = - |0)(0|. 

As follows from Eq. ((S]), the microwave pulse rotates 
the qubit state such that the initial density matrix pq 
becomes pq = URPqllj^, where 

Ur = e-'«"-^/2 = cos(6l/2) - in ■ asm{e/2) 

cos(6i/2) e-*^sin(6'/2) " 
-e^'^sin(6'/2) cos(6l/2) 

Here 9 ~ j^^ dtn{t), where ti and t2 are the pulse start- 
ing and ending time moments. 

The probability of the qubit to be found in state \i) is 

=Tr(|z)(z|p,) = Tr(P,p,) (z-0,1). (10) 

Here Pi is the projection operator Pi = ul^\i){i\Uji, i.e., 

p... 1 / l-fcos0 e-'-^sinf^X 
^i(")=2l e^-^sin^ 1-COS0 ) = 2^I + ^-'^)^ 

1 



(9) 



Po(a) = /-Pi(a)--(/-a-a), 



(11) 
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where / is the identity matrix and a is the unit vector 
a — (cos sin 9, sin </> sin 0, cos 9) (12) 

defining the measurement axis. 

Equations (fTO|) and (jlip show explicitly the equiva- 
lence between the qubit and detector rotations. Namely, 
one can interpret pi (po) as the probability of the qubit 
to be found in the state with the pseudospin parallel (an- 
tiparallel) to the measurement axis a. 

Notice that the microwave phase (f) is naturally defined 
modulo 27r, while the Rabi rotation angle 9 can be always 
reduced to a 2n range (we will assume — tt < 9 < it). 
Nevertheless, with this restriction there are still two sets 
of angles (0, (j)) corresponding to the same measurement 
direction a, because a is invariant under the change 

(61,0) <^ (-0,0 + tt). (13) 

It is easy to make a one-to-one correspondence between 
the measurement direction a and angles {9, 4>) by limiting 
either 6* or to a 7r-range (instead of 27r). However, we 
prefer not to do that because it is convenient and natural 
physically to have a 2tt range for one angle when the other 
angle is fixed. So, as follows from Eqs. (fT2|) and (|13p. the 
polar (zenith) and azimuth spherical coordinates of a are 
equal to, respectively, 9 and when 9 > 0, and —0 and 
(p + n when 9 < 0. 

The joint probability of the two-qubit measurement 
can be written as 

p.,,ia,b)^TT[FnS)P^ib)pl (14) 

where p is the two-qubit density matrix, P"" = Pi ® /, 
and Pi: ^I®P,. 

III. MAXIMUM BI VIOLATION: IDEAL CASE 

The purpose of this paper is to analyze conditions 
needed to observe the BI violation in experiment. Since 
it is usually easier to observe an effect when it is maximal, 
we start the analysis with the situations where violation 
of the BI is maximal. 



A. Bell operator and Cirel'son's bounds 

Equations © and (HH) yield E{d, h) = Tr(^_Bp), where 

A = P'^{d)-PS{a)^a-aa (15) 

and similarly B = b- Ub- Here (Ja = $ ® I and = / (E) (?, 
the eigenvalues of A and B being ±1. Correspondingly, 
as follows from Eq. 

5 = Tr(ep), (16) 
where the Bell operator^" B is 

B^AB-AB' + A'B + A' B' (17) 



(here A' = a' ■ Ua and B' = h' ■ db). 

The maximum and minimum values of S (so-called 
Cirel'son's bounds^) 

S± = ±2V2 (18) 

can be obtained, for example, in the following way4^ 
Since the Bell operator B is Hermitian, S± are equal 
to the maximum and minimum eigenvalues of B. These 
eigenvalues can be found analyzing the eigenvalues of B^: 

B^ = A+[A,A'][B,B'] = 4-4(axa'-CTa) (bxb'-db) (19) 

(here the vector product is taken before the scalar prod- 
uct). Since the eigenvalues of the Pauli matrices are equal 
to ±1, the largest eigenvalue of B^ is 8, achieved when 
a J- a' and & _L 6'. In this case the maximum and mini- 
mum eigenvalues of B are ±V8, thus leading to Eq. (|18p . 
(Both values ±-\/8 are realized because in this case the 
other eigenvalue of is 0, and the sum of all four eigen- 
values of B should be equal to since TtB ~ 0.) 

Consider some useful properties of S. As follows from 
Eqs. (|16p and (|17p . the value of S is invariant under ar- 
bitrary local unitary transformations Ua and Ub, 

p^{Ua®Ub)p{Ul®Ul), (20a) 

if simultaneously A UaAUl, A' -> UaA'Ul, B 
UbBUl, B' ^ UbB'ul or, equivalently, 

a RaS, a Raa , b Rbb, b' Rbb\ (20b) 

where Ra (Rb) is the rotation matrix corresponding to 
Ua (Ub), so that, e.g., Ua{a ■ d)U^ = {RaS) ■ ff. This 
invariance is an obvious consequence of the equivalence 
between the qubit and detector rotations, discussed in 
the previous section. As a result of the invariance, if some 
state is known to violate the BI for a given configuration 
of the detectors, one can obtain many other states and 
the corresponding detector configurations providing the 
same BI violation, by using Eqs. ([^0]) with all possible 
local rotations. 

Note also that B inverts the sign if the pair of vectors 
a, a' (or &, 6') inverts the sign. Correspondingly, for a 
given state 

S^-S if a ^ -a, a' -a' (or b -b, b' ^ -6'). 

(21) 

As follows from Eq. (|2ip . there is a one-to-one corre- 
spondence between the classes of detector configurations 
maximizing and minimizing S for a given state. 

B. Optimal detector configurations for maximally 
entangled states 

For any given detector configuration satisfying the con- 
dition 

d±a' and b±b', (22) 
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each of the Cirel'son's bounds (jlSp is achieved for a 
unique maximally entangled statei^i^ In contrast, for 
a given maximally entangled state there can be many 
optimal detector configurations giving the maximal BI 
violation. To the best of our knowledge, only the config- 
urations with the detector axes lying in one plane have 
been usually considered in the literature, though gener- 
ally the detector axes for different qubits may lie in two 
different planes. Moreover, the BI violation has been 
studied mainly for one of the Bell states;^ 

|*±) = (|10)±|01))/^/2, (23) 
|$±) = (|00>±|11»/V2, (24) 

while an arbitrary maximally entangled state can be writ- 
ten as 

l*me) = (|x?X?) + lx^X^))/%/2, (25) 

where {|xi): 1X2 )} orthonormal basis for qubit k. 

Our purpose here is to determine all optimal detector 
configurations for any maximally entangled state. 

1. Singlet state 

Let us start, assuming that the qubits are in the singlet 
state l^*-). For this state E{a,b) = (^"-1(0 • ct) (g) (6 • 
ct)|*_) = -a-b, so that [see Eq. ©] 

S = a-{b' -b)-a' -{b + b'). (26) 

Maximizing this formula over a and a', we should choose 
a to be parallel to b' — 6, while a' should be antiparallel 
to b' + b. Therefore"*^ a _L a', since b' — b and b + b' 
are mutually orthogonal. Similarly, rewriting S as S — 
b' ■ (a + a') — b ■ (a — a'), we can show that maximization 
of S requires & _L 6'. In this way we easily show that the 
necessary and sujjicient condition for reaching the upper 
bound S+ = 2^/2 for the singlet state is 

a_La', 6 -(a + a')/V2, 6' = (a' - a)/\/2. (27) 

Because of the symmetry (pij) . the necessary and suffi- 
cient condition for reaching the lower bound S- — —2^/2 
can be obtained by the inversion of the detector direc- 
tions for one of the qubits: 

a±a, b= {a + a)/V2, b' ^ {a- a)/V2. (28) 

Equations ((27| and ((28| show that for the singlet state 
the maximum BI violation requires that the detector axes 
for both qubits lie in the same plane. However, the orien- 
tation of this plane is arbitrary, since Eqs. ([TT]) and ([^5]) 
determine only the angles between the detector axes. 

All configurations maximizing (or minimizing) S for 
the singlet state can be obtained from one maximiz- 
ing (minimizing) configuration by all possible rotations 



of the plane containing the detector axes. As the ini- 
tial maximizing case we can choose the most standard 
configuration^i^ when all detector axes are within xz 
plane, 

<Pa - <P'a =<Pb= = 0, (29) 

and the polar (zenith) angles of the detector directions 
a, a', 6, b' are 

Oa^O, e'^^n/2, eb^-3TT/4, e[^-7r/4. (30) 

Then all detector configurations with S = 2\f2 can be 
parametrized by three Euler angles*® ki , K2 , and K3 (0 < 
^^1,3 "Si Stt, < K2 < "'), which describe an arbitrary 
rotation of the configuration P5)) ~ (P0)) . 

Similarly, all minimizing configurations [S — — 2\/2) 
can be obtained from the standard xz case [Eq. (|29p ] 
with 

9a = 0, 9'a^n/2, 9b = Tr/A, 9', = 37:/A (31) 

by arbitrary rotations of this configuration, characterized 
by three Euler angles ^1,2,3- 



2. General maximally entangled state 

Any maximally entangled two-qubit state can be ob- 
tained from the singlet state by a unitary transformation 
of the basis of one of the qubits^ (i.e. a one-qubit rota- 
tion). Therefore, an arbitrary case corresponding to the 
bounds S± = ±2^2 can be reduced to the singlet state 
considered above by a unitary transformation Ub of the 
qubit b basis and simultaneous corresponding rotation of 
the detector axes b and b' for the second qubit. Since 
the transformation Uf, can also be characterized by three 
Euler angles and K3, an arbitrary situation with 

s = 2V2 can be characterized by six independent pa- 
rameters {ki, K2, K3, k\, K2, K^), using the standard con- 
figuration ((29)) - (l30)) as a starting point. Similarly, any 
situation with S — —2\p2 is characterized by the same 
six parameters, starting with the xz configuration pip . 

Since these six parameters can describe arbitrary direc- 
tions of four measurement axes (a, a', 6, &') still satisfying 
the conditions a _L a' and 6 _L 6', it is obvious that any 
such four-axes configuration produces S = 2^2 for ex- 
actly one entangled state and also produces S = —2\p2 
for another entangled state. Notice that the sign of S 
can obviously be flipped by a 7r-rotation of qubit a (or 6) 
around the axis a x a! (b x b') instead of the 7r-rotation 
((2T|) of its detector axes. Also notice that six indepen- 
dent parameters for an optimal configuration can be al- 
ternatively chosen as any parameters characterizing the 
four measurement axes, which are pair-wise orthogonal: 
a _L a' and 6 _L 6' . 
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3. Odd states 

An important special case is the class of "odd" maxi- 
mally entangled states 

I*) = (|10)+e*"|01))/\/2 (0<a<27r), (32) 

which is of relevance for experiments with Joscphson 
phase qubitsji^ Such states can be obtained (with an 
accuracy up to an overall phase factor) from the sin- 
glet state 1 5*-) by unequal rotations of the two qubits 
around the z axis. Indeed, since Uzif) — e~'"^'^'^^ ro- 
tates a spin 1/2 around the z axis by angle tp, we obtain 
[U,{ao) O U^iao + tt - a)]|*-) = -ie-^"/^!*), where 
ao is arbitrary. Thus, in view of Eq. (PO]) . the optimal 
detector configurations for the odd state ([5^ can be ob- 
tained from those for l^"-) by rotating the detectors for 
the qubit b around the z axis by the angle n — a [notice 
that the state l|32p reduces to the singlet for a — n]. In 
terms of the parameters and 0, this is equivalent to the 
change 

'/>6 — > '/>6 + TT — a, (f)[ ^ (f)[ + TT — a. (33) 

Thus, for the odd states the class of optimal configura- 
tions maximizing S (as well as the class minimizing S) is 
characterized by four parameters: ki, K2, k^, and a. 

Now let us focus on the optimal configurations with the 
detector axes lying either in a "vertical" plane for each 
qubit (i.e., a plane containing the z axis) or the "horizon- 
tal" (xy) plane; such configurations will be important in 
the study of effects of errors (Sec. IIVP and decoherence 
(Sec. El. 

To obtain all "vertical" cases with the maximum BI 
violation S — 2\/2, we start with the standard configura- 
tion (p9)) - p0|) for the singlet, then apply a rotation in the 
xz plane by an arbitrary angle C (we can also apply the 
mirror reflection), then rotate the resulting configuration 
around the z axis by an arbitrary angle (po and finally 
apply the a-rotation ([55]) determined by the phase of the 
odd state ([5^ . As a result, the optimal measurement 
directions for the qubits a and b generally lie in different 
vertical planes, 

'/'a = 0a = </'0, 0& = (/^b = (/"O + TT - a, (34) 

while the polar angles corresponding to = 2\/2 are 

{9a, 9'a, 9,, 9',) = ±(0, 7r/2, -37r/4, -7r/4) + C, (35) 

where 0o a-nd C are arbitrary angles, while a is deter- 
mined by the state ((32|) . Notice that the signs ± corre- 
spond to the possibility of the mirror reflection, which 
we did not have to consider in the previous subsections 
because it can be reproduced using 3D rotations, while 
it is a necessary extra transformation in the 2D case. 

Similarly, the minimum S = — 2\/2 for the odd state 
([32)1 is achieved for the vertical configurations within the 
planes given by Eq. (|34p for the polar angles 

{0a,9'^, 9^, 9',) = ±(0, Tr/2, 7r/4, 3^/4) + C. (36) 



Recall that we define both 9 and cj) modulo 27r, and 
therefore each measurement direction corresponds to two 
sets of {9,(j)) [see Eq. (fT3| ]. Consequently, the optimal 
configurations described by Eqs. can be also 

described in several equivalent forms by applying the 
transformation p3|) to some of the four measurement di- 
rections. 

As follows from Eq. ([34|) . the only odd states for which 
the optimal vertical configurations lie in the same plane 
are the Bell states l^*-) (corresponding to a = tt) and 
1^1'+) (corresponding to a = 0). The optimal vertical 
configurations for the singlet state l^"-) are given by 

(l^a ^ (t>'a = (f'b ^ (f>'b ^ (1)0 (37) 

and Eqs. ([55]) To describe the optimal vertical con- 
figurations for the state it is natural to apply the 
equivalence (jl3[) to the qubit b measurement directions, 
so that the angles cj) are still all equal as in Eq. ([57)1 , while 
the angles 9 are given by Eqs. (|35 p - ([55|) with flipped signs 
for the qubit b, i.e. 9i, —9), and 9'^^ — > —9'j^. 

Now let us consider the optimal detector configurations 
in the horizontal {xy) plane: 

9a = 9'a = 9b = 9i = tt/2. (38) 

All configurations for S — 2^/2 can be obtained from the 
standard configuration (|29p - (PDl) by rotating it into the 
xy plane (so that the angles 9 are essentially replaced 
by the angles (p), then applying an arbitrary rotation 
within xy plane and possibly the mirror reflection, and 
finally applying the transformation (|33p with the state- 
dependent parameter a, so that 

{(t)a,(t)'a,<t>b + a, -f a) = ±(0, Ti/2, 7r/4, 3^4) + C (39) 

with arbitrary C (the signs ± correspond again to the 
possibility of the mirror reflection). 

Similarly, all horizontal conflgurations corresponding 
to S* = —2^/2 for the odd states are described by the 
angles 

(<^a, ^b a, 0; a) = ±(0, 7r/2, -37r/4, -7r/4) + C. 

(40) 

Notice that the application of the equivalence (fT3|) to 
all four measurement directions changes 7r/2 into — 7r/2 
in Eq. ([25]), while Eqs. and (gO]) do not change, since 
the corresponding 7r-shift of angles can be absorbed by 
the arbitrary parameter C. 



IV. LOCAL MEASUREMENT ERRORS 

In this section we consider the effects of local (inde- 
pendent) measurement errors on the BI violation. 
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A. Error model 

The probabilities of the measurement results for a sin- 
gle qubit can be written in the form 

1 

P^' = ^-'P™ = Tr(Q,p,), (41) 

r?j.=0 

where Pm are the probabilities which would be obtained 
by ideal measurements, Fim is the probability to find the 
qubit in the state \i) when it is actually in the state |to), 
and operator Qi — FmPq + FnPi contains the projector 
operators Pq.i [see Eq. ((TU)) ]. The operators Qi satisfy the 
same condition as the POVM measurement operators;^ 
namely, Qi are positive and Qq + Qi = 1. The condition 
Po^ + Pi^ — 1 implies that Fom + ^^im — 1- Hence, the 
matrix F has two independent parameters, which can 
be chosen as the measurement fidelities Fq = Fqq and 
Fi = Fii for the states |0) and so that 

Pa' = FoPno + (1 - Fi)pn, pf = (1 - Fo)poo + F^pn 

(42) 

(here pij are the components of the one-qubit density 
matrix after the tomographic rotation and < Fq i < 1). 
It can be always assumed that Fq + Fi > 1, since in 
the opposite case the measurement results can be simply 
renamed: ^ 1; as a consequence, max(i^o,^i) > 1/2- 

Using the assumption that measurement errors for each 
qubit can be considered independently of the errors for 
the other qubit, the measured probabilities for a qubit 
pair can be written in the form^ 

1 

Pf^^ E ^™^nP™« = Tr(g,"g5p), (43) 

■m,n— 

where F^^ is the matrix Fim for qubit k and Q^ — 
F^qP^ + F^\P^ [see Eq. dH])]. 

In this section we will discuss the condition for the BI 
violation as a function of measurement fidelities Fq and 
Fi- Sometimes we will limit the analysis to the case of 
equal measurement fidelities for both qubits, 

Fr = F^^F,, (44) 

however, the case of different measurement fidelities for 
the two qubits is also of interest. (Different fidelities 
are especially of interest when the qubits have different 
physical implementations. For instance, in the case of an 
atom-photon qubit pair— the detection efficiency for the 
atom is nearly 100%, whereas the photon-detector effi- 
ciency is significantly less than 100%.) Notice that for the 
Josephson phase qubits the trade-off between the fideli- 
ties Fq and F^ can be controlled in the experimen t^^i"^^ 
for each qubit individually by changing the measurement 
pulse strength. 

Several special cases of our error model have been 
previously discussed in the literature, starting with the 
CHSH paper.** For example, in the problem of the detec- 
tor loophole^^^ the CH inequality with Fq = Fq = 1 is 



often considered; then F^ is called the detector efficiency; 
both the cases = F^ (Refs. MM) and F^ ^ (Ref. 
Isoh have been considered. Let us also mention the effect 
of nonidealities on the BI violation considered for the ex- 
periments on two-photon interference.^!^ The situations 
of Refs. I27II28I formally correspond to the special case of 
our model with 

pa ^ pa ^ p^^ pb ^ pb ^ p^^ ^45) 

Then the product (2Fq — l)(2Ff, — 1) equals either the 
visibility^'' or the product of the visibility and the square 
of the signal acceptance probabilityi^ 



B. General relations for S 

The Bell operator (|17p can be generalized to the case 
of measurement errors. Inserting Eq. (|43)l into Eq. ([3]) 
yields E{a,b) — Tr{A]3p), where 

A =Ql-Q^ = F^- Fq" + (F(? + Ff - 1) a • (46a) 
B =F^-F^ + {F^ + F^-l)b-<Jb. (46b) 

Therefore S can be expressed as 

S = TiiBp) (47) 

via the modified Bell operator 

B ^ AB - AB' + A'B + A'B', (48) 

where A' and B' are obtained from A and B by replac- 
ing a and b with a' and b', respectively. Notice that B 
is a Hermitian operator and therefore in some cases it is 
useful to think about the measurement of S* as a measure- 
ment of a physical quantity corresponding to the operator 
B (even though this analogy works only for averages). 

It is rather trivial to show^ that in the presence 
of local measurement errors the Cirel'son's inequality 
15*1 < 2\/2 remains valid (this fact can be proven-* for 
any POVM- type measurement). Moreover, a stricter in- 
equality for jS'l [see Eq. ([5(7|) below] can be obtained, us- 
ing the method similar to that of Ref. [5l|. We will prove 
this inequality for all pure two-qubit states, p = 
which automatically means that it is also valid for any 
mixed state p. Using notation (O) = Tr(0/9) = {tp\0\ij) 
for any operator O, we start with the obvious relation 
1^1 =j{A{B - B')) + {A'{B + B'))\ < \{A{B - ^'))| + 
\{A'{B + B'))\. The next step is to apply the general 

inequality |(0i02)|^ < (OiOl)(0|02) to both terms in 
the sum (this inequality is the direct consequence of the 
Cauchy-Schwartz inequality \{ipi\tp2)\'^ < (V'i|V'i)(V'2|V'2) 
for the vectors jV'i) = and \tfj2) = O2 [■'/'))• In this 

way we obtain 

1^1 < ,J(A^){{B - B'y) + ^{A'^){{B + B'f) (49) 
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(notice that operators A, A', B, and B' are Hermi- 
tian). As the next step, we notice that the eigenval- 
ues of A (as weh as eigenvalues of A') are 2Fi — 1 
and 1 — 2Fq, which follows from Eq. and the 

fact that the eigenvalues of a • CTq are ±1. Therefore, 
(i^) < {2F-,^ - 1)2 and (A'^) < (2F-,, - 1)2, where 
i^^ax = inax(Fo'', F;f ); and so from Eq. (05]) we ob- 



tain |5| < (2i^,?_ - 1) 



{{B - B' 



{{B + B'Y) 



Next, since ^/xi + y/x2 < \/2{xi + X2) for any positive 
numbers xi and 2:2, and using the relation {{B — B')"^) + 
{{B + B'f) ^ 2(i?2 + we obtain the inequality 

l^l < 2(2^;?^^^ - l) \/(i?2 +B'2). Finally, using the re- 
lations {B^) < {2F^^^^ - 1)2 and (B'2) < (2^^,^^,, - 1)2, 
derived in a similar way as above, we obtain the upper 
bound 



|5| <2V2(2^^„"_-l)(2f^^,,-l) 



(50) 



This upper bound is generally not exact and can be 
reached only in the case when the errors are symmetric in 
both qubits [Eq. (gS])], leading to Eq. dSS]) below. While 
the bound (|50p depends only on the largest measurement 
fidelity for each qubit, our numerical results show that 
the exact bounds S± shrink monotonously with the de- 
crease of all fidelities, if the errors are small enough to 
allow the BI violation (see below). 

A useful expression for S can be obtained from Eqs. 
PS)) - (|T7|) by separating the terms for the ideal case: 

S = 2ed + 2^^da' • 4 + 2ee+^- + e+^+So, (51) 

where = + Fj^ - 1 , ^i'' - ^o*" > '^'o is the value of 

S in the absence of errors [Eq. pB]) ]. and Sk is the Bloch 
vector characterizing the reduced density matrix for the 
qubit k, i.e. pk = Trfc/^fcp — {I + Sk ■ (j)/2. Notice that 
Sa = Sb = for a maximally entangled state and there- 
fore the second and third terms in Eq. (|5ip may increase 
\S\ for nonmaximally entangled optimal states. That is 
why in the presence of errors the states maximizing and 
minimizing S are usually nonmaximally entangled3^ (see 
below) . 

Notice that in the presence of errors, S still preserves 
the invariance with respect to the local transformations of 
qubits and simultaneous rotation of measurement direc- 
tions described by Eqs. (|20p . The symmetry described by 
Eq. (PT|) (sign flip of S for the reversal of one-qubit mea- 
surement directions) is no longer valid; however, it can 
be easily modified by adding simultaneous interchange 
Fq ^ Fi of onc-qubit fidelities: 



S —S if a ■ 
S^-S if b- 



-a, a 
b' 



-a', FS ^ ; (52a) 
F^ ^ F^ (52b) 



Obviously, S does not change if both transformations 
(I52p are made simultaneously. As a consequence, the 
maximum S-^- and minimum (optimized over the qubit 



states and over measurement directions) are invariant 
with respect to simultaneous interchange of measurement 
fidelities 



pa ^ pa JTb 



(53) 



while the extremum values change as 5+ —S-, S- — > 
—S+ if only one-qubit fidelity interchange {Fq Ff or 
Fq ^ Fi) is made. (The corresponding optimal states 
obviously do not change.) 

In the presence of measurement errors the magnitudes 
of the maximum and minimum of S generally differ, 
5"+ \S-\- However, as follows from the latter sym- 
metry, 5'+ = \S-\ (as in the ideal case) if the two mea- 
surement fidelities are symmetric (equal) at least for one 
qubit: 



Fg" = Ff or Fg^ = F^ 



(54) 



If the fidelities are symmetric for both qubits [the sit- 
uation described by Eq. (|^5)) ]. then the expression for S 
given by Eq. (|?T|) becomes simple: S = (2Fa — l)(2Fh — 
l)So and directly related to the value Sq without mea- 
surement errors. Then the extremum values 



5± =±2\/2(2F, -1)(2F;,-1) 



(55) 



are obviously achieved for any maximally entangled state 
under the same conditions as in Sec. IIIII Correspond- 
ingly, the requirement for the fidelities for a violation of 
the BI is21^ 



(2Fa - 1)(2F6 - 1) > 2-1/2 _ Q 



(56) 



Notice that when the measurement fidelities for the both 
qubits are the same. Fa = Fb = F, Eq. ((56|l reduces to 
the threshold fidelity 



F > 0.5-1-2-^/'* « 0.920, 



(57) 



while if the measurement for one of the qubits is ideal 
(for example, Fq — 1), then the BI violation requires 



Ffc > 0.5 2-3/2 _ Q^g54_ 



(58) 



C. Analytical results for maximally entangled 
states 

Let us first analyze the extremum values of S for the 
class of maximally entangled states ([25)1 . Since in this 
case Sa — Sb — 0, we obtain from Eq. (j5ip that for maxi- 
mally entangled states 

S - 2(Ff - Fo") (Fi^ - Fo^) + (Ff + Fg - 1){F^ + F^-l)So 

... 

is directly related to the corresponding quantity Sq in the 

absence of errors. Therefore the extremum values of S 
for maximally entangles states are 



^± = 2(F{'-F°)(Fi''-Fo*)±2V2(Ff +Fo"-l)(Fi^+Fo^-l), 

(60) 
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and they are achieved under the same conditions as dis- 
cussed in Sec. IIIII 

When the asymmetry of measurement fidelities is sim- 
ilar for both qubits (i^f > and F^ > F^^ or both 
inequalities with "<" sign), the first term in Eq. ((60|) is 
positive, and therefore the BI jS*! < 2 can be stronger 
violated for positive S than for negative S. Similarly, if 
the asymmetries are opposite (for example, Ff > Fq and 
Fi < Fq), then it is easier to violate the BI for negative 
S. If fidelities are symmetric at least for one qubit [Eq. 
([51)1 ] . then the first term in Eq. (pO)) vanishes, and there- 
fore S- = —S+, as discussed in the previous subsection. 

If the fidelities are the same for both qubits, Fq — 
FI^ = Fo and Ff = F^ = Fi, then the positive S is 
preferable and 

S+ = 2(Fi - Fo)2 + 2x/2(Fi + Fo - 1)^ (61) 

(see the dashed lines in Fig. [T|) . This value of 5+ reaches 
the maximum 2\/2 when both errors vanish (Fo = Fi = 
1) and decreases with the decrease of each fidelity in the 
interesting region 5*+ > 2 (more accurately, as long as 
5'+ > 4-2v^ w 1.17). The condition for the BI violation 
in this case is 



(Fi - Fo)2 + V2{Fi +Fo- if > 1. 



(62) 



This threshold of the BI violation on the Fq-Fi plane is 
shown by the lowest dashed line in Fig. [1] It is an arc of 
the ellipse (corresponding to 5*+ = 2), which is symmetric 
with respect to the line Fq = Fi and is centered at Fo — 
Fi = 0.5. However, as seen from Fig. [1] this threshold 
looks quite close to a straight line on the Fq-Fi plane. 
Notice that in the case Fq = 1 the threshold reduces 
to the most well-known conditionii^ 



Fi > 2V2~2 « 0.828, 



(63) 



while in the case of symmetric error, Fq = Fi = F, we 
recover Eq. ((57|). 



D. Numerical results 

To optimize the CHSH inequality violation in presence 
of the measurement errors over all two-qubit states, in- 
cluding non-maximally entangled states, we have used 
numerical calculations. The analysis has been performed 
in two different ways (with coinciding results). First, we 
have searched for the maximum violation by finding ex- 
trema of the eigenvalues of the modified Bell operator B 
defined by Eq. (gS]). Since S is a Hermitian operator and 
S = Tr{Bp), for fixed measurement directions (a, a', b, b') 
the maximum and minimum eigenvalues of B are equal 
to the maximum and minimum values of S, optimized 
over the two-qubit states. Therefore, optimization of the 
eigenvalues over the measurement directions gives the ex- 
trema of S. Similar method has been previously used^^ 
for the case of identical local errors (H^ with Fq = 1, 




0.65 



0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 
Fo 

FIG. 1: Solid lines: contour plot of S+ (the maximum quan- 
tum value of S) versus the measurement fidelities Fo and Fi 
(assumed equal for both qubits), optimized over all two-qubit 
states. The dashed lines show the result of S+ maximization 
over the maximally-entangled states only [Eq. (|61|l ]. The Bell 
(CHSH) inequality can be violated when S+ > 2. 



while we apply this method to our more general error 
model. 

The numerical maximization (minimization) of the 
largest (smallest) eigenvalue of the Bell operator B has 
been performed using the software package Mathemat- 
ica. The full optimization should be over all four mea- 
surement directions (a, a', b, b'), described by 8 angles to- 
tal. However, because of the invariance of S under local 
transformations, it is sufficient to optimize B over only 
two angles: the angle between a and a' and the angle 
between 6 and 5', while the other angles are kept fixed. 

Our numerical results show that in general the optimal 
values of these two angles are different from each other; 
however, for equal fidelity matrices [Eq. ([44]) ] these angles 
are equal, so that a ■ a' = b ■ b' . This result has been 
obtained previously^^ for the special case Fq = 1. Our 
numerical results also show that 5*+ > |S'_| for positive 
values of the product (Ff - Fo'')(Fi^ - Fq*-) and 5'+ < |S'_ | 
when this product is negative, similar to the result for 
the maximally entangled states [see discussion after Eq. 

mi 

We have checked that the numerical results for S+ and 
\S-\ obtained via optimization of the eigenvalues of the 
Bell operator B coincide with the results (see Fig. [T|) ob- 
tained by our second numerical method based on the di- 
rect optimization of S. The second method happened 
to be more efficient numerically; as another advantage, 
it provides the optimal measurement directions together 
with optimal values 5+ and while the Bell-operator 
method gives only S+ and S- . 

In principle, direct optimization of S (for fixed mea- 
surement fidelities) implies optimization over the two- 
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qubit density matrix and over 8 measurement directions. 
However, there is a simplification. It is obviously suf- 
ficient to consider only pure states, since probabilistic 
mixtures of pure states cannot extend the range of S. 
Moreover, it is sufficient to consider only the states of 
the form 

I*) = cos(/3/2)|10) + sin(/3/2)|01), (64) 

since any pure two-qubit states can be reduced to this 
form by local rotations of the qubits (which are equiv- 
alent to rotations of the measurement directions); this 
fact is a direct consequence of the Schmidt decomposition 
theoremi^ The angle /3 can be limited within the range 
< P < TT because the coefficients of the Schmidt de- 
composition are non-negative. This range can be further 
reduced to < /3 < tt/2 since tt- rotation of both qubits 
about a;-axis (or any horizontal axis) exchanges states 
1 10) |01) and therefore corresponds to the transforma- 
tion /3 — > TT — /3. 

Our numerical optimization of S within the class of 
two-qubit states (|64p has shown that for non-zero mea- 
surement errors (we considered 2/3 < i^^*^ < 1) the opti- 
mal measurement directions (a, a', 6, b') always lie in the 
same vertical plane [such configuration is described by 
Eq. l|57|) ]. This vertical plane can be rotated by an ar- 
bitrary angle about z-axis [such rotation is equivalent 
to an overall phase factor in Eq. therefore we can 

assume </)o = in Eq. (|37|) . Notice that for the state 
([M)) the vectors Sa and Sf, in Eq. ([5T|) are along the z- 
axis, Sa = — Sb = i*cos/3. These vectors are zero for the 
maximally entangled state (/9 = 7r/2), then the vertical 
configuration is no longer preferential; however, the max- 
imally entangled state is optimal only when there are no 
measurement errors. 

For the state ([64)1 and vertical configuration ((37| of the 
detector axes the expression for S has the form 

-|-2cos/3(C+dcos6i;, -CC+cos6ifc), (65) 

where 

g = cos 9a cos 9b — cos 6a cos 9'ij + cos 9'^ cos 9i, 

+ COS 9'^ cos 9'ij, 
h = sin 9a sin 9f) — sin 9a sin + sin 9'^ sin 9f) 

+ sin 9'^ sin 91 (66) 

The numerical maximization and minimization of S in 
this case involves optimization over 5 parameters: /3, 6a, 
9'^, 9b, and 6'^. Nevertheless, in our calculations this pro- 
cedure happened to be faster than optimization over only 
two parameters in the method based on the Bell operator 
eigenvalues (we used Mathematica in both methods). 

The solid lines in Fig. [1] show the contour plot of max- 
imum value 5+ on the plane Fq-Fi for the case when 
the measurement fidelities for two qubits are equal [Eq. 
in this case S+ > \S-\]. Notice that the line for 



S+ = 2 ends at the points^^ Fq ^ 1, Fi = 2/3 and 
Fq = 2/3, Fi = 1 (strictly speaking, this line corresponds 
to 5*+ = 2 + since S — 2 can be easily realized with- 
out entanglement). The dashed lines, which correspond 
to the optimization over the maximally entangled states 
only [Eq. coincide with the solid lines at the points 

_Fo = -Fi, because in this case the optimum is achieved 
at the maximally entangled states, as follows from the 
discussion after Eq. ([54)l . When Fq ^ Fi, the use of 
non-maximally entangled states gives a wider range of 
measurement fidelities allowing the BI violation. How- 
ever, as seen from Fig.[Tl the difference between the solid 
and dashed lines significantly shrinks with the increase of 
S-i^ , so that there is practically no benefit of using non- 
maximally entangled states for the BI violation stronger 
than S+ > 2.4. Notice that the solid and dashed lines in 
Fig. [T] are symmetric about the line Fq — Fi since the in- 
terchange Fq ^ Fi does not change S+ , as was discussed 
after Eq. ([5^. 

Numerical calculations show that in the case of equal 
fidelity matrices, Eq. ([331), the optimal detector configu- 
rations for a nonmaximally entangled state (j64p have a 
"tilted-X" shape: a — —b' and a' = ~b. In this case the 
number of parameters to be optimized in (|65p reduced 
from 5 to 3, significantly speeding up the numerical pro- 
cedure. 

Notice that each optimal configuration within the class 
of states (l64l) corresponds to a 6-dimensional manifold 
of optimal configurations, obtained by simultaneous lo- 
cal rotations of the measurement axes and the two-qubit 
state (see discussion in Sec. IIIIB 2p . 



V. DECOHERENCE 

The detailed analysis of the effects of decoherence will 
be presented elsewhere In this sections we discuss only 
some results of this analysis, and also discuss the com- 
bined effect of local measurement errors and decoherence. 

To study effects of decoherence we assume for simplic- 
ity that the qubit rotations are infinitely fast. Thus, we 
assume that after a fast preparation of a two-qubit state p 
there is a decoherence during time t resulting in the state 
p' , which is followed by fast measurement of p' (including 
tomographic rotations). Now S is given by Eq. (jl6p (in 
the absence of errors) or (|47)) and ((5T|) (in the presence 
of errors) where p should be substituted by p' . To obtain 
p' we assume independent (local) decoherence of each 
qubit due to zero-temperature environment, described by 
the parameters 7fc = exp(— t/Tf ) and = exp(— </T2'^) 
(here k — a,b) where and ^^'^ the usual relaxation 
times for the qubit k {T^ < 2Tf ). 

As the initial state we still assume the state of the 
form ([Ml (even though in presence of decoherence this 
state actually does not always provide^ the extrema of 
S). It can be shown analytically^ that in the absence 
of measurement errors the maximum violation of the BI 
for the state ([M)) can be achieved when the detector 
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axes lie in either a horizontal [Eq. ([38)1 ] or vertical [Eq. 
([57)1 ] plane (any other detector configuration cannot give 
stronger violation). In the case of only population relax- 
ation {T2 = 2T^) the horizontal configuration is better, 
while in the case of only T2-effect {T^ — 00) the vertical 
configuration is better. 

When local measurement errors are considered to- 
gether with decoherence, the optimal detector configu- 
rations may be neither vertical nor horizontal. To eluci- 
date this fact, note that in Eq. ([^ with p replaced by p' 
the vectors Sa and S}, remain vertical in the presence of 
decoherence. As a result, when in the absence of errors 
the optimal configuration is horizontal, measurement er- 
rors may make the optimal detector axes to go out of the 
horizontal plane. Note, however, that for some parame- 
ter ranges the vertical and horizontal configurations are 
still optimal. 

In numerical calculations we should optimize S over 8 
parameters: /3 and 7 detector angles (one of the angles (p 
can be fixed because of the invariance of S under identi- 
cal rotations of the qubits around the z axis). We have 
performed such optimization for a few hundred param- 
eter points, choosing the measurement fidelities and 
decoherence parameters 7^ and Afc randomly from the 
range (0.7, 1). For many (more than half of) parame- 
ter points the optimal configuration was still found to be 
either vertical or (in much smaller number of cases) hori- 
zontal. Even when the optimal configuration was neither 
vertical nor horizontal, we found that restricting opti- 
mization to only the vertical and horizontal configura- 
tions gives a very good approximation of the extrema S± 
(within 0.01 for all calculated parameter points). Such 
restriction significantly speeds up the calculations, since 
we need to optimize over only 5 parameters instead of 8 
parameters. 

Assuming initial state ([64]) and replacing p by p' in Eq. 
(|T7)) we obtain 

S = + - 7a - 76 - (7a - lb) cosp]g 

+XaXbh sin P} + 2^^^ {ja + 7a cos /3 - 1) cos 0'^ 
+2ee+(7b-7bCos/3-l)cos0b, (67) 

when the detector axes are in a vertical plane [g and h 
are defined in Eq. ([M)) ]. and 

S = 2Cd + e+^l^h sin/3[cos(</)a - 0fc) 

- COs(0a - <l)'b) + COs((/)'„ - (jjb) + cos(0^ - 0^)] (68) 

for a horizontal detector configuration. 

To find the extrema of S within the class of vertical 
configurations, Eq. (p7|) should be numerically optimized 
over the parameter P and four angles 9. The optimiza- 
tion of S within the class of horizontal configurations is 
much simpler, because the term in the square brackets 
in Eq. can be optimized independently of /3. This 
optimization is exactly the same as in the ideal case [see 
Eqs. I|39p and ((40|) with a = 0], therefore the term in the 
square brackets has extrema ±2\/2, and therefore the 



extrema of Eq. (|68|) are reached at /3 = 7r/2 (maximally 
entangled state), thus yielding a rather simple formula 

S± = 2C-t ± 2V2e+C'+KXb, (69) 

which depends only on the r2-relaxation and measure- 
ment fidelities. 
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FIG. 2: Contour plot of the maximum value S+ versus the 
measurement fidelities Fq and Fi for the initial state (|64p in 
presence of decoherence with 7a = 7b = 0.96 and Aa = A;, = 
0.94, in absence (solid lines) or presence (dashed lines) of the 
symmetric crosstalk with pc — 0.1. Solid and dashed lines 
coincide for S+ — 2. 

For the numerical example shown by solid lines in Fig. 
[2] we assume that decoherence is identical for the two 
qubits and choose 7a = 7b = 0.96 and A^, = A;, = 0.94, 
that corresponds to realistically good values for phase 
qubits: Ti ~ 450 ns, T2 ~ 300 ns, and t ~ 20 ns. We also 
assume identical errors for both qubits, Eq. (j44p . which 
implies S+ > \S-\ (as in the absence of decoherence), so 
in Fig. [2] we show the contour plot only for S+. We have 
found that for these decoherence parameters the vertical 
detector configuration is better than any other configu- 
ration [assuming initial state Ij64p ] for any measurement 
fidelities in the analyzed range (0.8 < ^o,i < !)■ Notice 
that for the assumed decoherence parameters S+ = 2.50 
in the absence of measurement errors, and the BI viola- 
tion requires F > 0.947 (for Fq = Fi ^ F), that should 
be compared to the threshold F > 0.920 in absence of 
decoherence. 

Let us mention that the error model previously dis- 
cussed in relation to the BI violation in the two-photon 
interference^'^ can be shown to be formally equivalent 
to the special case of our model with pure dephasing 
{Ti =00) and identical errors with Fq = 1. Then our 
quantities Fi and AaA;, correspond, respectively, to the 
detector efficiency and visibility in Ref. [53. The case of 
pure dephasing in the absence of errors has been also con- 
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sidered in connection with the BI violation in mesoscopic 
conductors >2i 



VI. MEASUREMENT CROSSTALK 

The nonidealities discussed above are common for 
many types of qubits. Now let us discuss a more specific 
type of error: the measurement crosstalk for Josephson 
phase qubits. The crosstalk error originates from the 
fixed (capacitive) coupling between the qubits, which is 
still on in the process of measurement. The mechanism 
of the crosstalk is the foUowingp^ If the measurement 
outcome for the phase qubit a is "1" , then this qubit is 
physically switched to a highly excited state (outside of 
the qubit Hilbert space), and its dissipative oscillating 
evolution after the switching affects the qubit b. As a 
result, the extra excitation of the qubit b may lead to 
its erroneous switching in the process of measurement, 
so that instead of the measurement outcome "1,0" we 
may get "1,1" with some probability Similarly, be- 
cause of the crosstalk from the qubit b to the qubit a, 
we may obtain the measurement result "1,1" instead of 
"0,1" with some probability p'^. The values of p° and 
significantly depend on the timing of the measurement 
pulses applied to the qubits If the qubit a is measured 
few nanoseconds earlier than the qubit 6, then p\\ 
if the qubit b is measured first, then p° <C p^- the case 
when the measurement pulses are practically simultane- 
ous, the crosstalk probability becomes significantly lower 
and pj? « p^. 

Let us model the crosstalk in the following simple way. 
Even though physically the crosstalk develops at the 
same time as the measurement process and its descrip- 
tion is quite nontrivial^^ we will assume (for simplicity) 
that the crosstalk effect happens after the "actual" mea- 
surement (characterized by measurement fidelities as in 
Sections llVl and IV)) . so that the only effect of the crosstalk 
is the change of the outcome "1,0" into "1,1" with prob- 
ability p° and the change of the outcome "0,1" into "1,1" 
with probability Moreover, we assume that the prob- 
abilities pI: and p\ do not depend on the measurement 
axes (a, &, etc.). 

Notice that the measurement crosstalk obviously vio- 
lates the fundamental assumption of locality, on which 
the BI is based (so, strictly speaking the BI approach is 
not applicable in such situation). In this section we dis- 
cuss the modification of the classical bound for tak- 
ing the crosstalk into account (this bound is now model- 
dependent, in contrast to the usual BI), and we also an- 
alyze the effect of the crosstalk on the quantum result 
for S±. Besides that, we discuss a simple modification of 
the experimental procedure, which eliminates the effect 
of crosstalk by using only the "negative result" outcomes. 



A. Modified Bell (CHSH) inequality 

First, let us briefly review the derivation of the CHSH 
inequality, presented in Ref. [3^ . In a local realistic theory 

S = j s{K)F{K)dK, (70) 

where F{K) is a distribution of the hidden variable A and 

s(A) = A(A, a)S(A, 6) - A(A, a)B{k, V) 

+A{K, d')B{A, b) + A(A, a')B{K, b'). (71) 

Here the measurement outcomes ^(A, a) and B{A, b) can 
take only the values ±1, depending on the hidden vari- 
able A and the detector orientations for the qubits a and 
b. (Notice that in Sections |TTT] and [IV] the notation A 
and B has been used for operators; now they are clas- 
sical quantities. Also notice that the outcome value — 1 
is associated with the result "0".) It is easy to check 
that s(A) = ±2 under the locality assumption: the result 
A{K, a) does not depend on the orientations 6 of the qubit 
6 and vice versa; similarly, -F'(A) does not depend on a 
and b. After integration (iTOl) . this leads to the CHSH 
inequality H]). 

In our model the measurement crosstalk cannot change 
the positive product of outcomes AB = 1; however, it 
changes AB = —1 into AB = 1 with the probabil- 
ity if ^ = -B = 1 or with probability Pc(A) if 
A = —B = ^l. (The locality assumption is obviously 
violated, since the value of A now depends not only on 
A and a but also on B and thus implicitly on 6; similarly 
for B.) Notice that we assume that the crosstalk proba- 
bilities p° and may in principle depend on A (this is a 
slight generalization of the more natural A-independent 
model for p° and p^). 

The random change from AB = — 1 to AB = 1 due 
to the crosstalk leads to the modification of the CHSH 
inequality ([T]). Here we mention only the main points 
of the derivation of the modified CHSH inequality; the 
details are in the Appendix. We start with fixing A and 
considering all possible changes of the quantity s(A) in 
Eq. ((7T|l due to the crosstalk, thus obtaining the new val- 
ues of s together with their probabilities for each of 16 
realizations of the vector C = {A, A' , B, B'). It is easy 
to see that due to the crosstalk s(A) can get the values 
±4, which are outside of the limits ±2. Averaging the 
value s(A) over the crosstalk scenarios and then maxi- 
mizing and minimizing the result over 16 realizations of 
the vector C, we get 

-2+4minK(A),p^(A)} < (s(A)) < 2+2|p^(A)-p^(A)|. 

(72) 

Finally averaging this result over A, we obtain the mod- 
ified CHSH inequality: 

- 2 + 4minK,p^} < ^ < 2 + 2^^ - (73) 

which is the main result of this subsection. 
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Let us consider two special cases. For a symmetric 
crosstalk, = p[l — Pc, the inequality ((75|) becomes 

- 2 + 4pc < S' < 2, (74) 

while for a fully asymmetric crosstalk, pj! = 0, it becomes 

-2<5<2 + 2p°. (75) 

(similarly for = 0). 

It is interesting to notice that the inequality in the sym- 
metric case is more restrictive than the BI ([T]) ("easier" 
negative bound and no change of the positive bound). 
This is actually quite expected because in the limiting 
case Pc = 1 the crosstalk makes all AB products equal 1 , 
so that S = 2 always, as also follows from Eq. ([74|) . In 
contrast, for the fully asymmetric crosstalk the inequal- 
ity ((75)) is less restrictive than ^ ("harder to violate" 
positive bound and no change of the negative bound). In 
the case of a finite crosstalk asymmetry, both the positive 
and negative bounds change [Eq. ([73| ]. 

Let us emphasize that in contrast to the derivation of 
the original CHSH inequality, s(A) may be significantly 
outside of the range (—2, 2); it can get the values s(A) = 
±4 for any (symmetric or asymmetric) crosstalk. So the 
fact that the lower bound for S never decreases and the 
upper bound increases only slightly for small crosstalk, 
is due to the statistical averaging of the random increase 
and decrease of s(A) due to the crosstalk. 

B. Quantum calculation of S 

In the quantum case the crosstalk changes the mea- 
surement probabilities pij — > as 

Poo =Poo, Pio = (1 -Pc)Pio, Poi = (1 -Pc)Poi. 
=Pii+PcPio+Pcm- (76) 

Then using the definitions and ^ we obtain that the 
measured value of S becomes 

= 2pe + (1 - Pc)S + 2(p^ - pl)[pb{h) ~ Pa{d% (77) 

where Pc = (Pc +Pc)/2j while 5, Pa{a'), and Pb{b) are the 
quantities obtained in the absence of crosstalk [pa{a') and 
Pb{b) are defined after Eq. ([5])]. 

In a general case the maxima and minima of Eq. 
([77)1 can be found numerically. To estimate the effect of 
the crosstalk, let us calculate Eq. ((77)) for a maximally 
entangled state in the absence of errors and decoherence. 
Then pa{a') — pb{b) — 1/2, and using S± = ±2%/2 we 
obtain 

= 2V2 - {2V2 - 2)pc, ^ -2V2 + {2V2 + 2)pc. 

(78) 

As we see, both extrema are affected by the crosstalk, 
making the range narrower from both sides; however, the 
lower boundary is affected much stronger than the upper 



boundary. Comparing Eq. (|78)) with the modified CHSH 
inequality (|75)) . we see that the lower bound shifts up for 
the quantum result always faster than for the classical 
bound; therefore the gap between the quantum and clas- 
sical bounds always shrinks due to the crosstalk. The 
classical-quantum gap at positive S shrinks from both 
sides due to the crosstalk. 

In the case of symmetric crosstalk, p" — Pc — Pc, we 
can easily consider non-maximally entangled states, mea- 
surement errors and decoherence, since Eq. (|77)) in this 
case reduces to 5"" = 2pc -I- (1 — Pc)S'. Therefore, are 
simply related to the values S± without crosstalk (but 
with measurement errors and decoherence): 

S^ = 2pc + il-pc)S± (79) 

[similar dependence was used in Eq. (175)) ]. A violation of 
the upper bound of the modified CHSH inequality (|7il) 
can be observed when — 2pc -I- (1 — Pc)5'+ > 2, which 
yields 5*+ > 2, while a violation of the lower limit in 
Eq. dill) requires = 2pc + (1 - Pc)S- < -2 + 4pc, 
which yields S- < —2. Quite surprisingly, the symmet- 
ric crosstalk does not change the conditions for the BI 
violation. (Of course, the violation of the increased lower 
bound is not as convincing psychologically as the viola- 
tion on the increased upper bound.) 

Let us discuss the combined effect of local errors, deco- 
herence, and symmetric crosstalk for the numerical exam- 
ple considered in Sec. [V) assuming symmetric crosstalk 
with Pc = 0.1. Now for the state ([M)) in the absence of lo- 
cal measurement errors (Fq = Fi = 1) we get = 2.45, 
which is slightly less than the value S+ = 2.50 obtained 
in the absence of the crosstalk. The dependence of 
on the measurement fidelities Fo and Fi in this case is 
illustrated by the dashed lines in Fig. O In accordance 
with the above discussion, the solid line for Sj^ = 2 and 
the dashed line for — 2 coincide (the BI violation 
boundary is not affected), while the comparison of the 
solid and dashed lines for S+ = 2.2 and 2.4 shows that 
the crosstalk makes an observation of a given BI- violating 
value of S+ more difficult. 

Figure [3] shows a similar contour plot on the Fq-Fi 
plane for the lowest quantum value assuming the 
same parameters as in Fig. [2| Comparing Figs. [2] and [3] 
we see that it is more difficult to violate the lower classical 
bound (even though it is now only —1.6) than the upper 
classical bound of 2. This is because we assumed the 
same measurement fidelities for both qubits, that gener- 
ally shifts the quantum result up [as in Eq. ()60p ] . Notice 
that in Fig. [3] the reference bound of —2 can still be vio- 
lated, though only for almost perfect fidehties. 

C. Elimination of the crosstalk effect 

Slight modification of the CHSH inequality ((T)) can 
make it insensitive to the crosstalk. The main idea is 
to use only experimental outcomes with the result "0", 
when a qubit does not switch, and therefore the crosstalk 
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FIG. 3: Contour plot of the quantum minimum S- versus the 
measurement fidelities Fo and Fi (same for both qubits) op- 
timized over the states (|64p . We assume symmetric crosstalk 
with pc = 0.1 and decoherence parameters 7a = 76 = 0.96 
and Aa = Ai, = 0.94 (same as for Fig. The classical bound 
is shifted from —2 to —1.6 by the crosstalk. 

does not occur. Such "negative-result" ( "null- result" ) ex- 
periments with the Josephson phase qubits are very in- 
teresting from the quantum point of view^i^ however, 
here we are interested only in the classical consequence 
(or rather absence of it) for a measurement with a null 
result. 

Instead of the inequality ([T]) let us use the equivalent 
inequality ([5]), and let us change the definition of T in 
Eq. ([5]) by interchanging the measurement outcomes "1" 
and "0" . Then by symmetry the inequality ^ is still 
valid, so we get the classical bounds 

- 1 < f < (80) 

for 

T = pm{a,b) ~ pm{a, b') +poo{a',b) +poo{a',b') 

-poia')-poib), (81) 

where the probability po(a') is for measuring the qubit a 
only (without measuring the qubit 6), while poib) is for 
measuring only the qubit b. 

With this simple modification, the CHSH inequality 
becomes insensitive to the mechanism of the measure- 
ment crosstalki2i^ considered in this section. Notice, 
however, that in performing experiment in this way it is 
still important to check the absence of a direct crosstalk 
(due to the measurement pulse itself). This can be done 
by applying a measurement pulse to the well-detuned 
qubit b (so that it cannot switch) and checking that this 
does not affect switching probabilities for the qubit a 
(and similarly for a interchanged with b) . 



VII. CONCLUSION 

In this paper we have considered the conditions for the 
violation of the Bell inequality in the CHSH formi for the 
entangled pairs of solid-state qubits, when instead of the 
rotation of optical polarizers (detectors) we have to rotate 
the states of two qubits before the measurement, which 
itself is always performed in the logical z-basis (|0),|1)) 
for each qubit. While most of our results are applicable to 
many types of qubits, we have focused on the experiments 
with the Josephson phase qubits.*^ We have analyzed the 
BI violation for the ideal case as well as in presence of 
various nonidealities, including local measurement errors, 
local decoherence, and measurement crosstalk. 

In the ideal case the maximum violation of the BI 
{S± = ±2\/2, while the classical bound is \S\ < 2) can be 
realized for any maximally entangled state. The optimal 
configuration of the measurement directions in this case 
can be realized with three degrees of freedom for each 
maximally- entangled state (the measurement direction 
in our terminology actually refers to the qubit rotation 
before the measurement). However, in presence of non- 
idealities there is typically less freedom in choosing the 
optimal configuration. For the "odd" two-qubit states 
involving superpositions of the states |01) and |10), 
we have focused on the "vertical" measurement configu- 
rations, for which all measurement directions (a, a', b, b') 
are within the same vertical plane of the Bloch sphere, 
and the "horizontal" configuration, for which the four 
measurement axes are within the x-y plane. 

The qubit measurement with finite local errors (charac- 
terized by the fidelities Fq and Fi for each qubit) shrinks 
the quantum range for S. We have found that for a max- 
imally entangled state the BI violation is still possible 
when the classical bounds ±2 are exceeded by the ex- 
trema of the quantum result given by Eq. (|60p . In partic- 
ular, when two qubits have the same fidelities, the viola- 
tion condition is given by Eq. (j62p and shown by the low- 
est dashed line in Fig. [TJ it can be crudely approximated 
by the condition (Fi -|- Fo)/2 > 0.92. A significantly 
softer violation condition can be obtained when allow- 
ing the two-qubit state to be non-maximally entangled;^— 
this condition is shown by the lowest solid line in Fig. [T] 
However, the trick of using a non-maximally-entangled 
state does not help much when we need a BI violation 
with a significant margin (not just barely); this can be 
seen by comparing the solid and dashed lines in Fig. [T] 
For non-maximally entangled odd two-qubit states ((64l) 
in presence of local measurement errors, the vertical mea- 
surement configuration is found to be preferable in com- 
parison with other configurations. 

Analyzing the effect of local decoherence of the qubits 
for the odd two-qubit states ([M)) . we have found that ei- 
ther vertical or horizontal configuration of the measure- 
ment directions is optimal, depending on the parame- 
ters. In particular, in the case of population (energy) 
relaxation in 2-basis, the horizontal configuration is opti- 
mal, while for pure dephasing the vertical configuration 
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is optimal. In presence of both the decoherence and lo- 
cal measurement errors, the optimal configuration can be 
neither horizontal nor vertical; however, restricting opti- 
mization to only these two classes of configurations gives 
a very good approximation of the extrema S±. Obvi- 
ously, both the decoherence and the measurement errors 
make the observation of the BI violation more difficult. 

We have also analyzed the effect of the measurement 
crosstalki^ which plays an important role in measure- 
ment of the capacitively-coupled phase qubits. Since the 
crosstalk is a mechanism of classical communication be- 
tween the qubits, strictly speaking the BI is inapplica- 
ble. However, for a particular model of the crosstalk it 
is possible to derive a modified CHSH inequality [see Eq. 
(|73|1 ]. In particular, we have found that the symmetric 
crosstalk does not change the upper classical bound but 
increases the lower classical bound. The crosstalk also 
affects the quantum bounds, which are given by Eq. (j78[) 
for the maximally entangled state in the otherwise ideal 
case with arbitrary crosstalk and by Eq. (l79|) for an ar- 
bitrary case but assuming a symmetric crosstalk. Quite 
unexpectedly, the symmetric crosstalk does not change 
the threshold condition for the observation of the BI vi- 
olation. However, the crosstalk always reduces the gap 
between the classical and quantum bounds and makes an 
observation of the BI violation with a finite margin more 
difficult. It is important to mention that the detrimental 
effect of the crosstalk can be eliminated by a slight change 
of the CHSH inequality (by using only negative-result 
outcomes), which makes it insensitive to the crosstalk 
[see Eqs. §0^ and (|8T|) ]. 

We have performed the numerical simulations with the 
parameters similar to the experimental values for the 
best present-day experiments with the Josephson phase 
qubits. -i^iSS, Our results (see Fig. [2]) show the possibility 
of the CHSH inequality violation with a significant mar- 
gin even without further improvement of the phase qubit 
technology. 
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APPENDIX: DERIVATION OF THE 
INEQUALITY (pfsl) 

To derive the inequality ((72|) . we fix A and use the 
abbreviated notation s ~ s(A), p° — p^i^), Pc = PcW^ 
A = A(A, a), A' = A{A,a') and similarly for B and 
B' . Let us introduce the vectors C = {A, A' , B, B') and 
A = {AB',AB,A'B,A'B'), so that s ^ A - (-1,1,1,1). 
The vector C can assume 16 values, whereas without the 
crosstalk A can take 8 values, since the number of pluses 



or minuses in A is even. Each pair of C and —C yields 
one value of A. Generally, the crosstalk effect differs for 
C and —C, except for the symmetric case, p'^ = — pc, 
when crosstalk depends only on A. We start the analysis 
with the symmetric case and then consider the general 
asymmetric case. 



1. Symmetric crosstalk 

It is easy to check that without crosstalk s = 2 or 
—2. To obtain the upper bound of the modified BI, 
we consider four values of A corresponding to s = 2. 
The crosstalk cannot change A = (1,1,1,1), hence we 
discuss three other values: (—1, —1, 1, 1), (—1, 1,-1, 1), 
and (—1, 1, 1, —1). For any of these vectors, the crosstalk 
makes s to take the values 0, 2, and 4 with the probabili- 
ties pc^c, ql+pl, and pc^c, respectively, where qc = l-pc- 
Indeed, the change of any —1 to 1 in ^ occurs with the 
probability PcQc, yielding s = for the change of the 
first —I in A and s = 4 for the change of any other — 1, 
whereas +Pc is the probability of no change or change 
to A — (1, 1, 1, 1), both cases yielding s = 2. Thus, 
though now s can achieve the maximal mathematically 
possible value 4, it is easy to see that the maximal value 
for the average of s over crosstalk is still (s)inax — 2. 

To obtain the lower limit of the inequality, we consider 
four values of A corresponding to s = — 2 in the absence 
of the crosstalk. Consider first A = (— 1, — 1, — 1, — 1). 
This vector can be changed by the crosstalk to any of 16 
possible combinations of four pluses and minuses, yield- 
ing the values s — —4,-2,0,2,4 with the probabilities 
PcQc, 9c + 3Pc9c: ^Pcql + iplqc, iplql + pi, plqc, re- 
spectively (note that the sum of the above probabilities 
equals 1). The above probabilities can be easily obtained 
if one takes into account that s = —4 results from the 
change of only the first —1 in s — —2 occurs when 
either A have not changed or the first and one of the 
last 3 components have changed, s — Q occurs when the 
crosstalk results in A with the first and two other com- 
ponents equal to 1 or —1, s = 2 occurs when two of the 
last three components or all the components of A change 
the sign, and s — A results from the changes of all the 
last three components in A. As a result, in the case when 
only A = (—1, —1, —1, —1) is realized, we obtain 

(s)=-2 + 4p,. (A.l) 

Finally, let us consider the values A = (1,-1,-1,1), 
(1, 1, -1,-1), and (1,-1, 1, —1). For any of these vectors 
the crosstalk makes s to take the values —2, 0, 2 with 
the probabilities 5^, 2pcqc, and p^, respectively. These 
probabilities follow from the fact that s = —2 when A 
is not changed, s — results from the change of only 
one component —1, and s = 2 results from the change of 
the both negative components. It is easy to check that 
for any of the three above vectors, we again obtain Eq. 
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Combining the results for the upper and lower bounds, 
we obtain the inequality —2 + Apc < (s) < 2, the average 
of which over A yields the modified BI for the symmetric 
crosstalk given by Eq. ([7^ . 



two vectors C and —C corresponding to each A. As a 
result, the probability of each value of s obtained above 
is replaced by two probabilities, which differ from each 
other by the change p° ^ [cf. Eq. (|A.2[) ] . 



2. Asymmetric crosstalk 

As mentioned above, when p° ^ p^, the crosstalk yields 
different results for C and —C. However, recalling that 
the change AB = —1^1 occurs with the probability 
p^ if A — 1 and _B = — 1 or with the probability p^ if 
A = —I and B ~ I, we can obtain a symmetry relation 
for the probability of a change of A with a given C to 
some A' as a function of p° and p'^ : 

Pa^A'{p''c,pI,C) = Pa^A'{pIp'c,-C). (A.2) 

To extend the results of Sec. [T] to the general case of 
asymmetric crosstalk, we should reconsider all possible 
changes of A discussed in Sec. [l]with the account of the 



It happens that for most combinations A the prob- 
ability obtained in Sec. [T] is replaced by two proba- 
bilities just by replacing pc with or with p^. One 
of two exceptions is the case A — (—1,1,-1,1). Then 
for C = (—1,1,-1,1) the quantity s — 2 changes to 
s = 0,2,4 with the probabilities q^p''^, q^q^ + PcPc: 
and Pc^c: respectively, where = 1 — p^. This yields 
{s) — 2 + 2{p'^ — Pc)- The other value of (s) [correspond- 
ing to C = (1,-1,1,-1)] is obtained from this formula 
by substitution p'^ ^ Pc, so that (s) 2 -f 2(p^ - p^). 
The other exception isy^=(l,— 1,1,— 1). In this case 
we obtain (s) = —2 + 2{p'^+p'^) (for both possible values 
of C). Finally, by choosing the worst cases for the lower 
and upper bounds for (s) and averaging the result over 
A, we obtain Eq. ([75)1 . 
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